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Abstract We present the UV spectroscopy and timing 
of three nearby pulsars (Vela, B0656+14 and Geminga) 
recently observed with the Space Telescope Imaging Spec- 
trograph. We also review the optical and X-ray proper- 
ties of these pulsars and establish their connection with 
the UV properties. We show that the multiwavelengths 
properties of neutron stars (NSs) vary significantly within 
the sample of middle-aged pulsars. Even larger differ- 
ences are found between the thermal components of Ge- 
minga and B0656-I-14 as compared to those of radio-quiet 
isolated NSs. These differences could be attributed to 
different properties of the NS surface layers. 
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1 Introduction 

Optical through X-rays radiation from a typical isolated 
neutron star (NS) is expected to exhibit two components: 
thermal radiation from the NS surface and non-thermal 
radiation from the NS magnetosphere. The spectrum of 
non-thermal radiation can be described by a power-law 
(PL) model with a spectral index — 1 <^ a ^ (F^ oc ly"). 
This radiation is commonly interpreted as synchrotron 
emission from relativistic electrons/positrons accelerated 
by the electric fields near the NS surface and from sec- 
ondary particles produced in the pair cascades. The mag- 
netospheric radiation is intrinsically anisotropic and, there- 

O. Kargaltsev 

Pennsylvania State University 525 Davey Lab, University 
Park, PA 16802 
E-mail: oyklOO@psu.edu 

G. Pavlov 

Pennsylvania State University 525 Davey Lab, University 

Park, PA 16802 

E-mail: pavlov@astro.psu.edu 



fore, it shows strong pulsations. The non-thermal emis- 
sion completely dominates the multiwavelength spectrum 
in very young pulsars (e.g.. Crab and B0540-69). 

As the magnetospheric emission becomes fainter with 
increasing pulsar age, thermal emission becomes detect- 
able at T = P/2P ^ 10 kyrs, being seen as a "thermal 
hump" on top of the flat non-thermal spectrum. This 
thermal radiation is emitted from the NS surface which 
can be non- uniformly heated (e.g., due to anisotropic 
heat conductivity of the crust and bombardment by rel- 
ativistic particles). Indeed, in several cases at least two 
thermal components with different temperatures are need- 
ed to fit the X-ray spectrum (e.g., Pavlov et al. 2002; 
Zavlin & Pavlov 2004a; De Luca et al. 2005). The cooler 
thermal soft (TS) component is commonly interpreted 
as emission from the bulk of the NS surface. The hotter 
thermal hard (TH) component is often attributed to the 
NS polar caps (PCs) which can be additionally heated by 
the magnetospheric particles. As expected, the TS tem- 
peratures generally decrease with pulsar age. However, 
they do not fall onto a monotonically decreasing cool- 
ing curve, which suggests that NSs may have different 
masses (e.g., Yakovlev & Pethick 2004). In X-rays, the 
TS component usually shows weaker and broader pulsa- 
tions than the TH component. The TS pulsations can be 
caused by non-uniformities of the NS surface tempera- 
ture and by anisotropy of local emissivity in the strong 
magnetic field. 

Spectra of middle-aged pulsars (10 ^ t ^ 500 kyr) 
are particulary interesting because they often exhibit all 
the three emission components. The nearby middle-aged 
pulsars. Vela (t « 11 kyrs), B0656+14 (hereafter B0656; 
T « 110 kyrs) and Geminga (r « 340 kyrs), are relatively 
bright and have been extensively studied in both the op- 
tical and X-rays. Given their different ages, the three 
pulsars provide a representative cut through the popu- 
lation of middle-aged pulsars, allowing one to study the 
evolution of the pulsar properties. 

The X-ray spectra of these three pulsars can be de- 
scribed by a three-component, TS+TH+PL, model, with 
the thermal components modeled as BBs. However, such 
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Table 1 Basic properties of three middle-aged pulsars 
and three RQINSs with measured spin-down parameters 
and/or parallaxes (Kaplan & van Kerkwijk 2005ab; see 
also van Kerkwijk & Kaplan, these proceedings). 
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Geminga 


~ 200 


340 


0.24 
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fits result in substantially different model parameters in 
each case. 

The youngest Vela pulsar shows the highest Tts ~ 
1.2 MK emitted from a region with i?TS ~ 3.7 km. This 
Rts is substantially smaller than a typical NS radius, 
which might indicate that a large fraction of the NS sur- 
face is too cold to be seen in X-rays^. A factor of ten 
older B0656 has a lower temperature, Tts ~ 0.7 MK, 
and a much larger emitting area, i?TS ~ 12 km, close 
to that of a typical NS. The oldest of the three pulsars, 
Geminga, is even colder, Tts ~ 0.5 MK, and has the 
largest emitting area (i?TS ~ 13 km). 

The PL component, with ax ~ —1, is the strongest 
in the Vela pulsar, Lx,pl ~ 3 x 10^^ erg s'^ (in 0.2-10 
keV), becoming weaker in B0656 (ax ~ —0.5; Lx.pl ~ 
1.4 X 10"^^ erg s~^), and much weaker in Geminga (ax ^ 
-0.6; Lx,PL « 2.2 x lO^o erg s'^). This PL component 
is typically interpreted as synchrotron emission produced 
by relativistic particles in the pulsar margentosphere. 

The strength of the TH component also appears to 
decrease with pulsar age. Being relatively strong in the 
Vela pulsar (Tth ~ 6.3x10^^ ergs~^), it becomes weaker 
in B0656 (Tth ~ 3.5 x 10'^^ erg s~^) and decreases dra- 
matically in Geminga (Tth ~ 4.4 x 10^^ erg s~^). The 
sizes of the TH emission regions are close to the conven- 
tional PC radii [rpc = (27ri?^s/cP)i/2; see Table 1] in 
Vela and B0656, while in Geminga, Rtu ~ 46 m, is sur- 
prisingly small when compared to its rpc — 440 m. This 
may suggest that the surface of Geminga is more uni- 
formly heated and the "TH component" simply mimics 
a phase-dependent PL component^. However, regardless 
of the exact nature of the TH and PL components, the 

^ Alternatively, it may mean that the BE description of the 
thermal spectrum is inadequate, and a more realistic atmo- 
sphere model should be used. For instance, the Vela pulsar 
spectrum fits equally well by a two-component, NS atmo- 
sphere plus PL, model, which gives T « 0.7 MK at i? = 13 
km (Pavlov et al. 2001; see also Kargaltsev 2004). 

^ From a phase-resolved analysis of the XMM data, Jack- 
son & Halpern (2005) conclude that the X-ray spectrum of 
Geminga can be alternatively described by a BB+PL model 
with a phase-dependent PL slope. 
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Fig. 1 Unabsorbed NIR through FUV spectra of the 
Geminga, B0656 and Vela pulsars (top to bottom). The 
PL+BB fits are shown by the solid blue lines. The con- 
tributions of the individual components are also shown 
(blue dashed line for PL, blue dash-dotted line for BB). 
The dashed red line shows the extrapolation of the TS 
component fitted to the X-ray spectrum (see also Fig. 3). 
For B0656 we have also plotted the fit to the NUV-f FUV 
spectrum at the pulse minima (black dotted line). 

soft X-ray (0.3—2 keV) spectra of these three pulsars are 
clearly dominated by the thermal emission. 

Observing the TS radiation over a wide range of fre- 
quencies is particulary important. X-ray observations can 
only probe the Wein tail of the TS spectrum emitted 
from the NS surface. This information is insufficient to 
reconstruct the overall shape of the spectrum which can 
deviate from the simple black-body due to the pres- 
ence of an atmosphere or nonuniformity of surface tem- 
perature. Chemical composition of the atmosphere can 
dramatically affect the X-ray spectrum (Romani 1987; 
Pavlov et al. 1995; Zavlin, Pavlov & Shibanov 1996; 
Zavlin & Pavlov 2002), with light element atmospheres 
leading to a large Wien excess. Any surface temperature 
inhomogeneities will also complicate the spectrum, with 
hotter regions being increasingly important at higher fre- 
quencies. For these reasons, comparison of thermal X- 
ray emission with optical-UV emission on the Rayleigh- 
Jeans (R- J) side of the thermal hump is particularly valu- 
able. The challenge here is that the non-thermal mag- 
netospheric emission becomes increasingly dominant at 
longer wavelengths. 

The observations with ground-based telescopes and 
the Hubble Space Telescope (HST) have provided multi- 
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Table 2 STIS observations 



Pulsar 


Date 


Instrument 


Exposure, s 


A\ ^, A 


^^cgs 


/p^ % 


Geminga 


2002 Feb 27 


NUV-MAMA /F25SRF2 


11367 


1800-3000 


1.7 X 10" 


-IB 


40 




2002 Feb 26 


FUV-MAMA/G140L 


10674 


1155-1702 


3.7 X 10" 


-15 


65 


B0656+14 


2001 Sep 1 


NUV-MAMA/PRISM 


6791 


1790-2950 


2.6 X 10" 


-15 


67 




2001 Nov 16 


NUV-MAMA/PRISM 


12761 


1790-2950 


2.6 X 10" 


-15 


89 




2004 Jan 20 


FUV-MAMA/G140L 


4950 


1153-1700 


4.2 X 10" 


-15 


64 


Vela 


2002 May 28 


NUV-MAMA /F25SRF2 


2895 


1800-3000 


7.2 X 10" 


-IB 


87 




2002 May 28 


FUV-MAMA/G140L 


3060 


1153-1701 


8.2 X 10" 


-IB 


73 



Note. — * — instrument -|-filter passband; ^ - observed flux in the corresponding passband (for the NUV observations of 
B0656, the flux has been measured from the two observations combined); - intrinsic (corrected for the background) pulsed 
fraction. 



band photometry from ncar-IR (NIR) to ncar-UV (NUV) 
and showed that the NIR-optical spectra of the three 
pulsars are non-thermal, with a hint of a R- J component 
seen in Geminga and B0656+14 at A < 3000 A. The 
observations in UBVRIHJ filters show that the optical 
spectrum of the Vela pulsar remains very flat from NIR 
to UV, with ao ~ -0.1 (Shibanov et al. 2003). For the 
other two pulsars, the NIR-optical spectra are steeper, 
ao « -0.4 (Koptsevich et al. 2001; Shibanov et al. 2006). 

For all the three pulsars optical pulse profile measure- 
ments have been carried out. The Vela and B0656, show 
strong non-sinusoidal pulsations (GouifTcs 1997; Kern 
et al. 2003). Optical pulsations of Geminga were only 
marginally detected in the B band (Shearer et al. 1998). 

Despite the extensive coverage in the optical and X- 
rays, very few NUV observations had been carried out 
until recently, and the FUV parts of the spectra have 
remained completely unexplored. To fill this gap, we have 
undertaken an observational campaign with the MAMA- 
NUV and MAMA-FUV detectors of the Space Telescope 
Imaging Spectrograph (STIS) aboard the HST. 



2 Observations with STIS in far- and neeir-UV 

We have observed all the there pulsars in the NUV (1800- 
3000 A) and FUV (1150-1700 A) passbands (see Ta- 
ble 2). The data reduction and analysis are described in 
detail by Kargaltsev et al. (2005), Romani et al. (2005), 
and Kargaltsev & Pavlov (2006). Here we summarize the 
main results of these observations. 



2.1 Geminga 

In the FUV passband we obtained a low-resolution spec- 
trum of Geminga using MAMA-FUV with the grating 
G140L. The observed flux in the 1155-1702 A range is 
jFpuv = (3.72 ±0.24) x IQ-^^ ergs s"^ cm'^, correspond- 
ing to the luminosity Lpuy = 47rrf^^FUV = (1.78 ± 



0.11) X lO^^^rf^QQ ergs s ^. Fitting the specitruni with 
the absorbed PL model gives apuv = 1-43 ± 0.53 for 
E{B -V) = 0.03. The PL slope is close to that of the 
R-J spectrum, F,y oc u'^ , suggesting that the observed ra- 
diation is dominated by thermal emission from the NS 
surface. To estimate the NS surface temperature, we fit 
the absorbed BB model to the spectrum. Since the FUV 
frequencies are in the R-J part of the spectrum, the tem- 
perature values are strongly correlated with the radius- 
to-distance ratio (approximately, T (x (P/R"^). For a typ- 
ical NS radius, i? = 13 km, and the assumed distance 
d = 200 pc, the inferred temperatures arc 0.31 ± 0.01 
and 0.41 ± 0.02 MK, for E{B - V) = 0.03 and 0.07, 
respectively. 

The Geminga pulsar was imaged with MAMA-NUV 
using the broad-band filter F25SRF2. The NUV flux de- 
pends on the assumed spectral slope and extinction. For 
a plausible E(B— V)=0.03 and Q!nuv = 1, the unab- 
sorbed NUV flux is J^nuv = 2.2 x 10~^^ ergs cm"^ s~^ 
in 1800-3000 A. 

We used the 125 /is time resolution of STIS for timing 
analysis of both NUV and FUV data. We found statisti- 
cally significant pulsations with frequencies lying within 
(-0.3,-1-0.7) /xHz around the / = 4, 217, 608.6953 /xHz 
frequency estimated from multi-epoch timing observa- 
tions in 7-rays and X-rays (Jackson et al. 2002). The 
folded (source plus background) NUV light curve, plot- 
ted in Figure 2, shows one broad (FWHM « 0.8 of the 
period), flat-top peak per period, centered at the phase 
1.0. The most notable feature of the pulse profile is 
the narrow dip at ^ w 0.45. The pulsed fraction, defined 
as the ratio of the number of counts above the minimum 
level to the total number of counts in the light curve, 
is about 28%, which corresponds to the intrinsic source 
pulsed fraction /p w 40%. 

The FUV light curve, plotted in the same Figure 2, 
also shows a sharp, asymmetric dip at (p « 0.45 with an 
additional shallower dip at ^ « 0.95. The pulsed fraction 
in the observed (source -|- background) radiation is about 
40% (intrinsic pulsed fraction 60%-70%.). 
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2.2 B0656+14 

The FUV spectrum of B0656 has been measured with 
MAMA-FUV using the G140L grating. Unfortunately, 
only two of the eight planned orbits of data were col- 
lected because of the STIS failure on August 3, 2004. The 
observed flux in the 1153-1700 A band is ^fuv = (4.3 ± 
0.3) X 10""'^^ ergs cm~^, corresponding to -Lpuv = 
(4.2 ± 0.3) X lO^^rf^gg ergs s^^ The fit with the ab- 
sorbed PL model results in a rather steep spectral slope 
apuv = 1-51 ± 0.62 for a plausible E{B - V) = 0.03. 
Similar to Geminga, this suggests that the FUV radia- 
tion could be dominated by thermal emission from the 
NS surface. The absorbed BB fit gives the surface tem- 
perature of 0.71 ± 0.03 MK for E{B -V) = 0.03, i? = 13 
km, and the distance of 288 pc (Brisken et al. 2003). 
The corresponding unabsorbed bolometric luminosity is 
Lho\ = (3.1 ± 0.5) X 10^^ ergs s"^ Thus, the measured 
brightness temperature is substantially higher than that 
of Geminga. 

The NUV spectrum and light curve of B0656 have 
been obtained using MAMA-NUV with PRISM (Shibanov 
et al. 2005). We re-analyzed these data to facilitate the 
direct comparison with the FUV and other multiwave- 
length data. The observed flux in the 1790-2950 A range 
is .Fnuv = (2.63 ± 0.38) x 10"^^ erg s"^ cm'^, corre- 
sponding to Lnuv = (2.61 ± 0.38) x lO'^^d'^^^ erg s^^. 
Fitting the spectrum with the absorbed PL model, we 
flnd onuv = 1.09±0.41 for a plausible E{B-V) = 0.03. 
The PL slope is somewhat flattc^r than that of the FUV 
spectrum, whicli can be explained by a larger contribu- 
tion of the non- thermal component. 

Making use of the event time-tags, we created FUV 
and NUV pulse profiles folded with the radio ephemeris 
provided by M. Kramer (2005; priv. comm.). The FUV 
pulse profile (Fig. 2) shows two peaks of approximately 
equal strengths. The minimum at « 0.55 is shallower 
than the other minimum at (p ~ 1.05. The observed 
pulsed fraction, « 36%, corresponds to the intrinsic puls- 
ed fraction of w 64%. We also folded the light curves 
for each of the NUV visits and co-added them into the 
combined NUV pulse profile (Fig. 2). The NUV pulse 
shape resembles the FUV one but shows a higher pulsed 
fraction (the observed and intrinsic pulsed fractions are 
22% and 80%, respectively). The lower pulsed fraction 
in FUV is evidence of a larger contribution of thermal 
emission at the FUV wavelengths (thermal emission is 
intrinsically less pulsed than the non- thermal emission). 
The NUV and FUV light curves can be interpreted as 
two non-thermal peaks on top of a single low-amplitude 
thermal pulse with the maximum at </) « 0.5. This inter- 
pretation would imply a substantially softer spectrum at 
the pulse minima than at the peaks. 

To verify this assumption, we produced the phase- 
resolved NUV through FUV spectrum for three phase 
intervals [0.5; 0.7], [0.9; 1.2], and [0.2; 0.5] [0.7; 0.9], cor- 
responding to each of the two peaks and two minima 



combined. The difference in the slopes of the PL model 

fitted to the peaks spectra does not exceed the spectral 
index uncertainties. However, the spectrum at the pulse 
minima is substantially steeper {a = 1.5 ± 0.2) than the 
peak spectra (a = 0.6±0.2). This behavior is in line with 
the expected larger contribution of thermal emission at 
the minima phases. 



2.3 Vela pulsar 

The spectroscopic (FUV) and imaging (NUV) data has 
been acquired using exactly the same instrumental setup 
as in the Geminga observations (see Table 1) and reduced 
in a similar way. 

For a plausible PL slope, qnuv = (sec §3.3) and 
color index E{B — V) — 0.05 (estimated from the hy- 
drogen column density found from the X-ray fits; San- 
wal et al. 2002), the unabsorbed NUV flux is ^nuv ~ 
1.1 x 10"" ergs cm-2 s"^ in 1800-3000 A. 

The FUV spectrum of the Vela pulsar is significantly 
flatter than those of Geminga and B0656, which reflects 
its predominantly non-thermal nature even throughout 
the FUV passband. The total observed flux in the 1153- 
1701 A range is J^^fuv = (8.19 ± 0.36) x lO'^^ ergs s"! 
cm-2, corresponding to Xfuv = (8.78 ± 0.39) x 10^^ 
ergs s-i at d = 300 pc (Dodson et al. 2003). Fitting the 
spectrum with the absorbed power-law model, we found 
apuv = 0.06 ± 0.39 at E{B -V) = 0.05. 

Figure 2 shows the NUV and FUV light curves folded 
with appropriate radio ephemeris (courtesy of R. N. Man- 
chester). At least four narrow peaks can be identified 
in each of the light curves. The higher signal-to-noise 
(S/N) NUV light curve shows some substructure in the 
two main peaks. The observed (source plus background) 
pulsed fractions arc 52% and 73% for the FUV and NUV 
light curves, respectively. The corresponding intrinsic pulse 
fractions are 73% and 87%. 



3 Multiwavelength spectra and light curves 

To compare the UV emission with the NIR, optical and 
X-ray emission, we plotted multiwavelength spectra of 
the three pulsars in Figure 1 (NIR to FUV) and Fig- 
ure 3 (NIR to X-rays). The FUV, NUV and X-ray light 
curves are shown in Figure 2. Below we briefly discuss 
the multiwavelength properties for each of the pulsars. 



3.1 Geminga 

The top panel of Figure 1 shows the UV fluxes together 
with the multi-band NIR-optical photometry adopted 
from Kargaltsev et al. (2005) and Shibanov et al. (2006). 

As the NIR through FUV spectrum cannot be de- 
scribed by a simple PL model, we fit this spectrum with 
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Fig. 2 NUV, FUV and X-ray (top to bottom) light curves of the Gemigna, B0656-f 14 and Vela pulsars (left to 
right). 



a BB+PL model. For the fixed R/d = 13 km/200 pc, 
we obtain T = 0.30 ± 0.02 MK, ao = -0.46 ± 0.12, 
-Fo = 0.11 ± 0.02 ^Jy for E{B - V) ^ 0.03, where ao 
and Fq are the parameters of the PL component: F^, = 
Fo (v/l X 10^^ Hz)"o. Notice that the parameters of the 
BB component are virtually the same as obtained from 
the FUV-MAMA spectrum alone. We see that the BB 
emission dominates at A ^ 3000 A, while the magneto- 
spheric PL emission dominates at longer wavelengths. 

To investigate the connection between the optical/UV 
and X-ray properties, we use the data obtained with 
XMM-Newton. The X-ray spectrum can be reasonably 
well described by a three-component TS-f TH+PL model 
(see §1 and Kargaltsev et al. 2005 for details), shown in 
the top panel of Figure 3. One can see that the extrapo- 
lation of the TS component overpredicts the FUV fluxes 
by a factor of 1.6 [for E{B - V) = 0.03; see also Fig. 1]. 
This FUV deficit can also be demonstrated by plotting 
together the temperature-radius confidence contours for 
the TS component and for the BB fit to the FUV-MAMA 
spectrum (Fig. 4). We see that at plausible values of in- 
terstellar extinction, E{B — V) < 0.07, the FUV contours 
lie at smaller radii (or much lower temperatures) than 
the X-ray contours, in contrast to some other neutron 
stars (see §4.2). The extrapolation of the X-ray PL com- 
ponent suffers from large uncertainties and is marginally 
consistent with the optical fiuxes. 

The background-subtracted X-ray light curves are 
compared with the UV and optical light curves in Figure 



2. In the 0.2-0.5 and 2-8 kev bands, the radiation is dom- 
inated by the TS and PL components, respectively The 
PL hght curve (pulsed fraction /p = 34% ± 8%) shows 
two pronounced peaks per period, resembling the 7-ray 
light curve (Jackson et al. 2002), albeit with a smaller 
distance between the peaks, and a hint of a third peak, 
at (/) « 0.2. On the contrary, the 0.2-0.5 keV light curve 
(/p = 30% ± 2%) is characterized by one broad peak per 
period (with small "ripples" , perhaps due to contribution 
from the PL and TH components). The minimum of the 
TS light curve is approximately aligned in phase with 
one of the minima of the PL light curve, being shifted 
by Z\0 « 0.1 from the sharp dips of the NUV and FUV 
light curves. 



As the thermal soft emission dominates the spectrum 
from 4 eV to 1 keV, the no n- thermal pulsations can only 
be observed in the optical or hard X-rays. The double- 
peaked shape of the hard X-ray (2-8 keV) pulsations 
does not correlate with the TS pulsations. The optical 
pulsations, marginally detected in the B band (Shearer 
et al. 1998), also show a double-peaked structure, with 
a larger separation between the peaks [Acf) k, 0.6 in the 
optical vs. « 0.4 in the hard X-rays). The phases 
of the optical peaks differ from those of the hard X-ray 
peaks. To understand and resolve these apparent incon- 
sistencies, a higher quality optical light curve should be 
obtained. 
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14 15 16 17 18 

Log i>, Hz 

Fig. 3 Unabsorbed multiwavelength spectra of the Geminga, B0656+14 and Vela pulsars (top to bottom). In all 
cases the X-ray spectra (thick red lines) are fitted with the TS+TH+PL model. The TS and PL X-ray components 
are shown with dashed red lines (the uncertainties of the PL fit are dotted lines). The TH component is shown by 
the green dashed line. The solid blue lines show PL+BB fits to the optical through FUV spectra; the blue dashed 
and dash-dotted lines are the contributions of the PL and BB components, respectively. 



3.2 B0656+I4 

Similar to Geminga, we fit the NIR through FUV spec- 
trum of B0656 (middle panel in Fig. 1; NIR and optical 
data are from Koptsevich et al. 2001 and Shibanov et al. 
2006) with a two-component, BB+PL, model. For the 
fixed R/d = 13 km/288 pc, we found the temperature 
T = 0.50 ± 0.05 MK, spectral index ao = -0.41 ± 0.08, 
and PL normalization Fq = 0.27 ± 0.02 ^Jy for E{B - 
V) = 0.02. One can see that the BB emission dominates 
at A ^ 2000 A. The fit also confirms that, in the NUV, 
the fraction of thermal photons is smaller than the frac- 
tion of nonthermal photons by a factor of 1.5, which is 
compatible with the deeper minima and higher pulsed 
fraction in the NUV light curve. 

We plot the NIR to X-rays spectra for B0656 in the 
middle panel of Figure 3. The X-ray spectrum is again 
fitted with the TS-f TH-fPL model. We see that the ex- 
trapolation of the TS component from X-rays to lower 
frequencies roughly coincides with the FUV spectral flux 
and goes slightly above (by a factor of 1.3) the UV ther- 
mal component. This excess, however, is not statistically 
significant (see Fig. 4). 

The double-peaked structure seen in the optical (Kern 
et al. 2003) and NUV/FUV light curves changes over to 



a single peak in the X-ray light curves (Fig. 2). It is 
possible that the asymmetric 0.25—0.7 keV pulse pro- 
file consists of two components: a broad one, peaked at 
about (f) « 0.60-0.65, and a narrower component cen- 
tered at « 0.90-0.95, close to the radio pulse phase, 
= 0. In this case the maximum of the broader com- 
ponent would approximately coincide with the shallower 
minimum in the NUV and FUV light curves, which is 
located at the phase of a presumable maximum in the 
thermal UV component (see §2.2). The peak of the nar- 
row component is approximately aligned with the single 
peak of the hard (2—4 keV) X-ray light curve. 

Given the non-thermal interpretation of the peaks in 
the optical- UV light curves and the close match between 
the optical and X-ray PL components (Fig. 3), it is sur- 
prising that the 2—4 keV pulse profile of B0656 shows 
only one peak (at the phase of the second optical/UV 
peak; (j) ~ 0.85). Since the phase-resolved spectroscopy 
shows similar PL slopes in both UV peaks, they are likely 
to be of a similar origin. The absence of the second peak 
at the X-ray frequencies could be explained assuming 
that the X-rays are emitted closer to the NS in a more 
narrow cone while the UV emission is generated higher 
in the magnetosphere and spread over a broader beam. 
Alternatively, it could be that the "X-ray PL" is not 
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Fig. 4 Top: Confidence contours (68%, 90%, and 99%) in 
the (soft) tcmpcraturc-radius plane obtained from fitting 
the EPIC-pn spectra (solid lines) with the TS+TH+PL 
and BB+PL models (labels near the contours). The 
dashed and dash-dotted lines show the confidence con- 
tours obtained from fitting the FUV spectrum with a 
BB model for two values of E{B ~ V). See Kargaltsev et 
al. (2005) for other details. Bottom: Similar confidence 
contours but for B0656-I-14. 



the other hand, the optical- UV spectrum is substantially 
flatter than the non-thermal X-ray spectrum, ax ~ — 1. 

The UV spectrum of the Vela pulsar allows one to ob- 
tain a restrictive upper limit on the NS surface temper- 
ature (see Romani et al. 2005 for details). The spectrum 
remains largely nonthermal even at the pulse minima, 
which constrains the NS temperature to ^ 0.4(d/300 pc)^ 
(i?/13 km)~^ MK, assuming the BB model for the ther- 
mal emission and E{B — V") = 0.05. The limit becomes 
less restrictive, ^ 0.6((i/300 pc)2(i?/13 km)-^ MK, if 
the hydrogen NS atmosphere model is used to describe 
the thermal spectrum. The inferred upper limit is sur- 
prisingly low compared to the temperature of the older 
B0656 (Table 3). 

In the UV, the Vela pulse profile is dominated by 
multicomponent nonthermal emission. In Figure 2 we 
show the NUV, FUV and X-ray light curves together. 
The X-ray light curves obtained with Chandra (Sanwal 
et al. 2002) and RXTE (Harding et al. 2002) are ex- 
tremely complex, with multiple peaks likely originating 
from different regions of magnetosphere or corresponding 
to physically different spectral components. In agreement 
with Harding et al. (2002), we identify at least five com- 
ponents in the multiwavelength pulse profile. No peaks 
are seen in the UV at the phases </> = 0.12 and 0.56 of 
the two 7-ray peaks (Kanbach et al. 1994). Instead, in 
the FUV they are replaced by tighter spaced peaks 1 
and 2, which continue to the optical. We also find that 
the UV peaks 2 and 4 have counterparts in the soft X- 
ray light curve while peaks 1 and 3 appear to be shifted 
by A(j) « 0.1 with respect to the other two soft X-ray 
peaks. The UV light curves show strong correlation with 
the hard X-ray light curves obtained with RXTE above 
10 keV (Harding et al. 2002), i.e. peaks 2, 3 and 4 are 
almost at the same phases in both UV and RXTE bands. 
At the same time, the correlation between the UV and 
2.5-8 keV light curves is less obvious. The spectral in- 
dices of the individual UV components are rather uncer- 
tain and show no correlation with the RXTE components 
(Romani et al. 2005). Should a better NUV/FUV data 
be ever obtained, a joint phase-resolved analysis of the 
UV and hard X-ray data may be more revealing. 



magnetospheric emission (sec §4.3), and the match be- 
tween the X-ray and optical spectra is just a coincidence. 



3.3 Vela pulsar 

Figure 1 shows the spectrum of the Vela pulsar from NIR 
to UV, with the optical/NIR data points adopted from 
Shibanov et al. (2003). This very flat spectrum nicely 
fits a PL model. The best-fit spectral index and nor- 
malization are ao = 0.01 ± 0.02 and = 1.50 ± 0.03 
liJy (at f = 10^^ Hz), respectively. The spectral index 
is larger than that obtained by Shibanov et al. (2003), 
ao = —0.12, from fitting only IR and optical data. On 



4 Discussion 

4.1 Thermal emission from NS surface 

The spectra of all the three pulsars show thermal emis- 
sion, but the strength and position of the thermal hump 
vary. In the Vela pulsar, thermal emission is seen only 
in soft (0.3-2 keV) X-rays while in Geminga and B0656 
the FUV spectra are also predominantly thermal. The 
observed FUV and soft X-ray spectra represent the R-J 
and Wien tails of the TS component. The TS BB tem- 
peratures measured in X-rays decrease from « 1.2 MK 
in the Vela pulsar to » 0.7 MK in B0656 and » 0.5 
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MK in the oldest Geminga. At the same time, the bolo- 
metric TS luminosity of the Vela pulsar is surprisingly 
low, which can be formally explained by a smaller emit- 
ting area of the TS component. The actual reason for 
this difference could be different properties (particulary, 
chemical composition and degree of ionization) of the 
NS atmosphere. In both Geminga and B0656 the bright- 
ness temperatures derived from the PL+BB fit to the 
combined NIR through UV data are lower than the tem- 
peratures of the X-ray TS component (see Table 3). The 
observed thermal UV spectrum of Geminga lies a factor 
of 1.5— 2 below the continuation of the X-ray TS compo- 
nent (for a plausible extinction), whereas the UV deficit 
is substantially smaller, if present at all, in B0656. In 
the case of Vela pulsar, we could only estimate an up- 
per limit on the UV brightness temperature, which lies 
above the extrapolation of the X-ray TS component. 

Turning to the light curves (Fig. 2), one can see that 
in Geminga and B0656 thermal soft X-ray light curves 
show broad maxima with some non-thermal "bumps" on 
top of them. The bump (centered at « 0.9) is relatively 
large in B0656, making the soft X-ray light curve asym- 
metric. The soft X-ray pulse profile of Vela is much more 
complex, reflecting the larger non-thermal contribution 
than in the other two pulsars. However, one still can ar- 
gue that there is a single broad thermal maximum in the 
(j) = 0.7 — 1.2 range. 

The FUV light curves of B0656 and, especially, Gemin- 
ga are also expected to be largely thermal. Yet, Geminga 
shows strong double-peaked FUV pulsations, different 
from both the NUV and the soft X-ray pulsations (see 
Fig. 3). Neither FUV nor soft X-ray pulsations of Gemin- 
ga can be produced by locally isotropic blackbody emis- 
sion. Although the soft X-ray pulsations can be explained 
by a magnetized atmosphere model (e.g., Zavlin & Pavlov 
2002), the high pulsed fraction and the peculiar shape of 
the FUV pulsations are not predicted by this model (no- 
tice, however, that the current atmosphere models are 
not applicable at the optical- UV frequencies) . A possible 
explanation of the strong thermal pulsations invokes a 
"screen" of absorbing plasma suspended in the NS mag- 
netosphere, which may partially eclipse the surface emis- 
sion at certain rotation phases. 



4.2 UV emission of Geminga and B0656 versus RQINSs 

Over the last decade seven soft X-ray sources with ther- 
mal spectra have been discovered by ROSAT. Very large 
X-ray-to-optieal flux ratios/limits and a lack of radio 
emission allow one to conclude that these objects are 
radio-quiet isolated neutron stars (RQINSs). Extensive 
X-ray observations of RQINSs with Chandra and XMM- 
Newton have provided high-quality spectra and period 
measurements for most of them (F. Haberl, these pro- 
ceedings). Recent period derivative measurements by Ka- 
plan & van Kerkwijk (2005ab) in RX J0720.4-3125 and 



RX J1308.6-h2127 (J0720 and J1308 hereafter) place these 
RQINSs near the anomalous X-ray pulsars (AXPs) in the 
P-P diagram. However, the X-ray properties of RQINSs 
are quite different from those of AXPs (the BB temper- 
atures of RQINSs are much lower, and no PL tails are 
seen). On the other hand, the X-ray spectra of RQINSs 
are also very different from those of ordinary old radio 
pulsars (t ^ 1 Myr; B ^ 10^^ G; e.g., Zavlin & Pavlov 
2004b; Kargaltsev et al. 2006). In terms of their ther- 
mal X-ray properties, RQINSs resemble middle-aged pul- 
sars (see Table 3); however, the optical-UV properties of 
RQINSs show important differences. 

Out of the seven RQINSs, five have been detected in 
the optical, and UV spectra have been obtained for two of 
the five, RX J1856.5-3754 (J1856) and J0720 (van Kerk- 
wijk & Kulkarni 2001; Kaplan et al. 2003). The UV spec- 
tra of these RQINSs look thermal {F^ v"^), and even 
their optical fluxes seem to follow the same R-J curves, 
although Kaplan et al. (2003) report a faint non-thermal 
PL compontent, with ao ~ 0.3, ipL.o ^ 6 x 10^^ ergs 
s-i for J0720 (for the other three RQINSs, the spectral 
slopes in the optical are uncertain) . We note that no de- 
tectable optical magnetospheric component is expected 
at least in J0720 and J1308 (and likely in other RQINSs, 
whose spin-down powers, E, have not been measured 
yet) if the nonthermal optical efficiency of these objects 
is similar to those of radio pulsars, Lpi^^o/E <^ 10~^- 
10~^ (Zavlin & Pavlov 2004b). For instance, one would 
expect LpL.o ~ 10^^-10^^ ergs s'^ for J0720, i.e. a fac- 
tor of > 10'' higher efficiency is needed to explain the re- 
sult of Kaplan et al. (2003). More importantly, RQINSs 
exhibit significant excess in thermal UV emission, with 
the UV fluxes exceeding the continuations of the X-ray 
thermal spectra by a factor of 5-9. The nature of this 
UV excess is not firmly established yet. It could occur if 
the UV and soft X-ray thermal components are emitted 
from different regions (with different areas) on a non- 
uniformly heated NS surface (Pavlov et al. 2002) . Alter- 
natively, these NSs may have a condensed (solid) surface, 
possibly covered by a tenuous atmosphere that is opti- 
cally thin in X-rays but optically thick in the optical-UV 
(e.g., Motch et al. 2003). Whatever is the nature of the 
large UV excess in RQINSs, we do not see it in B0656 
and see a UV deficit in Geminga. 

The apparently smaller UV-emitting area of Geminga, 
as compared to the X-ray-emitting area, cannot be ex- 
plained by a nonuniform temperature distribution. We 
might speculate that the temperature distribution over 
the bulk of Geminga's surface is more uniform than in 
J1856 and J0720, e.g., because of different geometry and 
strength of the magnetic field that affects the surface 
temperature distribution. However, to explain why the 
more uniformly heated Geminga exhibits quite substan- 
tial pulsations of its thermal X-ray radiation while no 
pulsations have been detected from J 1856, one has to 
assume very special orientations of the J1856's spin and 
magnetic axes. The lack of UV excess in Geminga can 
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also be attributed to different chemical composition or 
lower temperature of the Geminga's surface (Kargalt- 
sev et al. 2005). However, the latter argument does not 
apply to B0656 whose UV and soft X-ray thermal com- 
ponents show higher temperatures, comparable to those 
of J1856 and other RQINSs. One could also speculate 
that the surface of RQINSs is solid while in B0656 and 
Geminga it is in a gaseous state. Indeed, the magnetic 
fields of Geminga and B0656 are 5-10 times smaller than 
those of J0720 and J1308, implying lower condensation 
temperatures (van Adelsberg et al. 2005). In addition, 
irradiation of the NS surface by energetic particles and 
photons generated in the magnctosphcrc may ablate the 
NS surface and facilitate formation of thick atmospheres 
in active pulsars while we see no indications of magne- 
tospheric activity in RQINSs. The presence of thick at- 
mospheres (lack of solid sufrace) in pulsars is supported 
by the fact that a hydrogen atmosphere model provides 
a reasonably good description of the TS spectrum in the 
Vela pulsar while the BE model gives a very small NS ra- 
dius (Pavlov et al. 2001). (Such fully ionized atmosphere 
models, however, are not directly applicable to colder 
Geminga and B0656 because their atmospheres are not 
fully ionized.) In RQINSs, the solid surface, directly seen 
in X-rays, would emit the spectrum that resembles a BE 
but with the emissivity reduced by a factor of ^ 2 (van 
Adelsberg et al. 2005), which would explain the appar- 
ently smaller X-ray radii of RQINSs {Rts in Table 3). 
However, even a tenuous atmosphere on top of the solid 
surface can be opaque in UV, resulting in a larger UV 
emitting area. Within this interpretation, one still needs 
to explain the absorbtion features observed in the spec- 
tra of several RQINSs. Although it was suggested that 
these features can be formed in a hydrogen or helium 
atmosphere (see van Kerkwijk & Kaplan, these proceed- 
ings), they could also be due to the resonances at various 
hybrid frequencies in the dense magnetized condensate 
(van Adelsberg et al. 2005). 



4.3 Nonthermal emission 

In Geminga and E0656, the PL components fitted to 
the hard X-ray and optical spectra show similar slopes 
ax ~ Q^o ~ —0.5. However, in B0656 the extrapola- 
tion of the X-ray PL matches the optical points much 
better that in Geminga (see Fig. 3). For the Vela pul- 
sar, the slope of the PL component is much steeper in 
X-rays, ax ~ —1, than in the UV-optical, ao ~ 0-0. 
However, one should keep in mind that the slopes of the 
X-ray PL components are measured in the narrow 2-8 
keV band (Fig. 3) and hence arc very uncertain because 
of a large background at these energies and, possibly, 
because the spectrum deviates from a simple PL. For 
instance, a comptonized tail of thermal radiation could 
mimic the nonthermal PL component in the narrow 2- 
8 keV band. Measuring hard X-ray spectra above ^ 10 



keV will help to better understand the nature of the non- 
thermal emission in Geminga and E0656. 

The 2-8 keV light curves show rather sharp and strong 
peaks for all the three pulsars. If the 2-8 keV emission 

is indeed duo to the comptonization of thermal emission, 
one needs to explain its high anisotropy. The number of 
hard X-ray peaks varies from a single peak in E0656's 
light curve to a double-peaked structure in Geminga, 
and even more complex multi-peak light curve in the 
Vela pulsar. 

In the UV-optical, the Vela pulsar shows the strongest 
pulsations, with at least four narrow peaks. In contrast 
to the Vela pulsar, the UV light curves of B0656 have a 
large thermal contribution (« 70% in FUV and « 50% in 
NUV). Its NUV light curve shows two large peaks with a 
much wider separation than the main peaks in the NUV 
pulse profile of the Vela pulsar. Due to a much lower S/N 
in the B0656 NUV light curve, smaller and narrower non- 
thermal peaks may remain undetected. Similar to the 
Vela pulsar, the NUV light curve of E0656 is noticeably 
different from its 2 4 kcV light curve. No obvious non- 
thermal contribution to the UV light curves is seen in 
Geminga, which is consistent with a large thermal frac- 
tion in the UV. However, non-thermal pulsations should 
be better seen in the optical (A <j 4300 A). 



5 Conclusions 

We have observed three pulsars (Geminga, E0656+14, 
and Vela) at the far-UV and near-UV wavelengths with 
the STIS MAMA detectors and measured the pulsar 
spectra and pulse profiles. We have also analyzed the 
X-ray data and found that optical through X-ray spec- 
tra consist of thermal and non- thermal components in all 
the three pulsars. In particular, we found the following. 

1. Thermal contribution to the FUV passband in- 
creases with pulsar age. The thermal component is not 
seen in the younger Vela pulsar (r w 11 kyrs), but it 
dominates the FUV spectrum in the Geminga pulsar 
(r « 340 kyrs). As the FUV spectrum grows more ther- 
mal, the "thermal hump" , whose maximum is located in 
soft X-rays/EUV, shifts toward lower frequencies. 

2. The X-ray spectra of the two older pulsars, Gemin- 
ga and E0656, resemble those of RQINSs. However, their 
UV properties are noticeably different. In sharp contrast 
to RQINSs, whose spectra show large UV/optical ex- 
cess, Geminga shows UV deficit while the continuation 
of the X-ray TS component approximately matches the 
R-J component in E0656. Likely, this reflects the differ- 
ences in the magnetic field and phase state of the NS 
surface layers. 

3. While no UV pulsations have been detected in 
RQINSs, Geminga shows strong, non-sinusoidal pulsa- 
tions in the FUV range, where the spectrum is domi- 
nated by the thermal component. To explain the strong 
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Table 3 X-ray, UV and optical properties of middle-aged pulsars and RQINSs. 



Name 


IlJV 


log Lfuv 


log ipL.O 


ao 


Tts 


Rts 


log Lts 


Ttii 


Rtu 


log Ltii 


log LpL,X 


ax 




A IK 


(Tgs/S 


crgs/'s 




A IK 


km 


crgs/s 


A IK 


km 


orgs/'s 


crg/s 




Vela 


^ 0.4 


29.15 


28.57 


-0.01 


1.17 


3.7 


32.26 


2.1 


0.7 


31.80 


31.49 


-1.0 


B0656+14 


~0.53 


28.74 


27.87 


-0.41 


0.71 


12.2 


32.43 


1.4 


1.1 


31.54 


31.15 


-0.5 


Geminga 


0.30 


28.33 


27.17 


-0.46 


0.49 


12.9 


31.83 


2.3 


0.05 


29.64 


30.34 


-0.6 


RX J0720 


0.37 


28.44 


27.78? 


0.3? 


1.0 


5.1 


32.58 












RX J1856 


0.80 


28.76 






0.72 


8.7 


32.65 













Note. — First column is the brightness temperature measured from the TS+PL fit to the NIR through FUV spectra 
(see Fig. 1) for a 13 km NS radius. Second column gives the unabsorbed FUV luminosity in the 1150-1700 A passband 
[E{B — V)= 0.05, 0.03 and 0.03 for Vela, B0656 and Geminga, respectively]. Third and forth columns are the PL component 
slope and luminosity (in 4000-9000 A), respectively. Columns three through ten list the parameters of the TS-I-TH-I-PL fits 
to the X-ray spectra. Thermal component luminosities arc the bolometric luminosities while the unabsorbed PL luminosity, 
LpL,x = 47rd^JFpL,x, is in the 0.2-10 keV band. The luminosities and radii are calculated for the distances listed in Table 1. 



FUV pulsations, one may need to invoke magnetospheric 
absorption at certain rotation phases. 
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